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The deterministic Landau-Lifshitz-Gilbert equation has been used to investigate the nonlinear 
dynamics of magnetization and the specific loss power in magnetic nanoparticles with uniaxial 
anisotropy driven by a rotating magnetic field, generalizing the results obtained for the isotropic case 
found in [P. F. de Chatel, I. Nandori, J. Hakl, S. Meszaros and K. Vad, J. Phys.: Condens. Matter 
21, 124202 (2009)]. As opposed to many applications of magnetization reversal in single-domain 
ferromagnetic particles where losses must be minimized, in this paper, we study the mechanisms of 
dissipation used in cancer therapy by hyperthermia which requires the enhancement of energy losses. 
We show that for circularly polarized field, the loss energy per cycle is decreased by the anisotropy 
compared to the isotropic case when only dynamical effects are taken into account. Thus, in this case, 
in the low frequency limit, a better heating efficiency can be achieved for isotropic nanoparticles. The 
possible role of thermal fluctuations is also discussed. Results obtained are compared to experimental 
data. 
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I. INTRODUCTION 

The nonlinear dynamics of the magnetization in single- 
domain ferromagnetic nanoparticle systems has been the 
subject of an intense study and it is at date a chal- 
lenging issue. Examples are ferromagnetic resonance, 
switching of magnetization, data storage based on mag- 
netic devices, spintronics etc. The applications which 
are strongly related to the present work are the follow- 
ing: ferrofluids, magnetic resonance imaging (MRI) and 
other biomedical applications, see e.g. [lKD]. While in 
most cases the loss energy per cycle has to be minimized, 
in cancer therapy by hyperthermia the goal is to enhance 
the heating efficiency of magnetic nanoparticles driven 
by an external magnetic field, preferably inside the ma- 
lignant tumours. The common practice is to use a lin- 
early polarized external magnetic field alternating at a 
frequency of the order of 10 5 Hz. Indeed, in case of the 
linearly polarized applied field, the optimization of loss 
energy with respect to the amplitude and frequency of 
the external field has been studied in detail [4|, |8|, |9( . It is 
natural to ask what is the dependence of the specific ab- 
sorption rate on the nature of polarization, i.e. whether 
a better heating efficiency can be achieved by a circu- 
larly polarized applied field 0-0, EH ■ The study of dy- 
namical effects of circularl y p olarized field has received a 
considerable attention fll| - [l5j . Power losses for isotropic 
nanoparticles under rotating field have also been investi- 
gated in the presence (see e.g. p| @) or in the absence 
(see e.g. [l(|) of thermal effects, but no systematic analy- 
sis have been performed in order to investigate the effect 
of anisotropy on the energy absorption of nanoparticles 
in the low frequency limit suitable for hyperthermia. 

The relaxation and the loss energy of a single isotropic 
magnetic nanoparticle has been considered under circu- 
larly polarized applied field in fi"o| when no thermal ef- 
fects were included. In the low frequency limit, the loss 



energy per cycle was found to be larger in case of the lin- 
early polarized applied field as compared to the circularly 
polarized one. Thermal effects for isotropic system were 
studied in detail in |5|. In the limit of low frequency, 
the linearly polarized field was found to produce more 
heat power for higher temperatures, too. However, re- 
cent experimental results [7| show a different picture; the 
linearly and the circularly polarized external field pro- 
duced an equal heat power at least for low frequencies. 
Since, the immobility and the aggregation of particles 
into chains is a known feature of ferrofluids when the 
sample becomes very anisotropic, it is a natural question 
to ask whether the anisotropy can be responsible for the 
discrepancy between the theory and the experiment. 

The goal of this paper is twofold. On the one hand, we 
consider the role of anisotropy in the possible enhance- 
ment of heat power of magnetic nanoparticles (in the 
absence of thermal effects) by generali zing the results ob- 
tained for the isotropic case found in (l0(. On the other 
hand, we study whether the anisotropy can be used to 
explain the experimental results of Q • The possible role 
of thermal fluctuations is also discussed. 

The paper is organized as follows. In Sec|Hl the deter- 
ministic Landau-Lifshitz-Gilbert equation has been given 
in case of uniaxial anisotropy suitable for the description 
of magnetization dynamics for single magnetic nanoparti- 
cles (in the temperature range far from the Curie temper- 
ature) . The specific loss power and loss energy is studied 
in case of the circularly polarized applied field in Sec Mil 
Known results of the isotropic case is briefly summarized 
and compared to the findings of the present work done 
for nanoparticles with uniaxial anisotropy. In Sec IIV1 we 
study the linearly polarized applied field in the limit of 
large anisotropy and in Sec|V]the possible modification of 
the findings by thermal fluctuations is discussed. Finally, 
Sec I VII stands for the summary. 
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II. LANDAU-LIFSHITZ-GILBERT EQUATION 

In order to study the energy losses under repeated mag- 
netization reversal, one can distinguish two different pro- 
cesses related to the mobility of magnetic particles in 
fcrrofluids. Either the magnetic moment rotates within 
the particle (Neel regime [ljl) or the particle rotates as 
a whole (Brown regime). Another way to classify various 
types of relaxation mechanisms is related to the temper- 
ature. For example, if one considers relaxation far from 
the Curie temperature then the magnetization process in 
a single-domain particle can be well described by means 
of the Landau-Lifshitz [l7j equation which is mathemat- 
ically equivalent to Gilbert's one [ 1 8| with the appropri- 
ate definition of its coefficient [19|, [2fJ . This is referred 
as the Landau-Lifshitz-Gilbert (LLG) equation. For a 
complete description of relaxations close the Curie tem- 
perature, thermal effects should be incorporated. How- 
ever, in some cases (such as a rotating applied field) the 
LLG equation provides us reliable results on energy losses 
even at higher temperatures. 

In this work the focus is on the relaxation based 
on dynamical effect obtained by the LLG equation for 
anisotropic nanoparticles. We argue that at least for 
small and relatively large anisotropy, findings of the 
present work can be used to study energy losses in a tem- 
perature range relevant for hyperthermia. An important 
feature of the LLG equation is that the magnetization 
vector's magnitude does not change under the influence 
of the external held. Thus, it is convenient to rewrite it 
in terms of the unit vector M = m/mj where ms is the 
saturation magnetization. Then the LLG equation reads 
as 

^■M = - 7 '[M x H off ] + a'[[M x H off ] x M]. (1) 
dt 

with the coefficients 7' = /io7o/(l + ct 2 ) and a' = 
/ i o7o 7 ? TO s/(l + a2 ) where /io is the permeability of free 
space, 70 is the gyromagnetic ratio and rj is the dimen- 
sionfull and a — HoJoV m S the dimensionless damping 
constant. Let us introduce an effective gyromagnetic ra- 
tio 7 = 70/(1 + a 2 ). Then the parameters of the LLG 
equation (fTJ) can be rewritten as 7' — /io7 and a' = ^ja. 
It is important to note that the effective gyromagnetic 
ratio 7 used in this paper is positive as opposed to the 
negative parameter of [lfj . The cross denotes the vector 
product and the effective magnetic field acting on the 
magnetization M is defined as 

H c ff = H cxt + H aniso (2) 

with the alternating or circulating applied (external) field 
H oxt and the anisotropy field H an i so . Let us note, that in 
this work all considerations have been done in the absence 
of a static field. 



III. CIRCULARLY POLARIZED APPLIED 
FIELD 

In this section we discuss the solution of the LLG 
equations (fT]) obtained for an immobile single-domain 
(isotropic and anisotropic) magnetic particle under cir- 
cularly polarized, i.e. rotating applied (external) mag- 
netic field. The applied field is assumed to rotate in the 
xy-plane with an angular frequency to 

H GXt = —7 (cos(wi),sin(wi),0), (3) 
7 

where uj^ — 7'|H| is the Larmor frequency. (The angular 
velocity vector is perpendicular the the xy plane.) For 
the sake of simplicity we consider particles with uniaxial 
anisotropy where the easy axis of the magnetization is 
chosen to be the z-axis, i.e. the anisotropy field is defined 
as 

H aniso = ^ (0,0,X cS M z ), (4) 

7 

where M z is the z-component of the magnetization vec- 
tor and the parameter A e fi describes the strength of the 
anisotropy. With this particular choice of the anisotropy 
field the arrangement used in the paper is identical to 
that investigated in [ll], [ijj • For a more detailed anal- 
ysis, specially when the effect of thermal fluctuations is 
also taken into account, the easy axis of magnetization 
has to be chosen arbitrarily (for linearly polarized case 
see e.g. [22j]) but this is out of the scope of the present 
work. 

It is convenient to use a coordinate system in which 
the steady state solution of the LLG equation (fTJ) has the 
simplest form: a time-independent magnetization vector 
[10h12| . The transformation is done by an appropriate 
rotation [l(| 

(+ cos(cji) + sin(wt) \ 
-sin(u;i) +cos(wi) (5) 
lj 

which transforms the LLG equation into a coordinate 
system, which rotates around the z axis with the applied 
magnetic field. The transformed z axis points then in the 
direction of the angular velocity vector w. Denoting the 
Cartesian coordinates of the transformed magnetization 

(u Xl u v ,u z ) = OiM, (6) 
the LLG equation (fTJ) can be written as 

— ^- — LOUy + a N u 2 + ctNU 2 z 

-uj L \ c gUyU z - a N \ cS u x ul, 

dUy 

—— = -UJU X - UJ L U Z - OtNU x U y 

dt 

+UJl^cSU x U z - CCNKsUylll, 

= uj L u y - a N u x u z + ajvAoff (1 - u 2 z )u z , (7) 
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where = a'|H| is introduced. Let us note that A e ff 
is dimensionless but apf, uj and u>l are of dimension 1/s. 
However, by introducing a dimensionless time t = t/to all 
the frequency parameters can be rewritten as dimension- 
less quantities such as ajy = ajv^o etc. For the sake of 
simplicity we keep the original notation (without the tilde 
superscript) but the time and consequently the frequen- 
cies are considered as dimensionless parameters. Since 
the LLG equation retains the magnitude of the magneti- 
zation vector, only two of the Cartesian components are 
independent (u is a unit vector in the rotating frame). 
In order to describe the orientation of the magnetization 
let us introduce angles following the definition of 



u x = sin cos < 
u y = — sin 9 sin ( 



cos 9. 



(8) 



Thus in the rotating frame the LLG equation obtained 
for the three Cartesian coordinates (J7]) reduces to a set 
of differential equations for the two angles 



d0 

dt 

d0 

dt 



u>l sin0 + q;jv cos0 cos</> — a^Xeg sin 9 cos 9, 

uj L X cS cos 9. (9) 



cos V sin 

Uj, COS (f>— — - + w — Cep? 



sm( 



suit 



Let us note that the differential equations for 9 and <f> m 
© are identical to Eqs.(2,3) of [11| in case of vanishing 
static field (h az = 0) if one makes the following identi- 
fications lol = h a ±/(l + a 2 ), ajv = ah a ±/(l + a 2 ) and 

A c ff = Kes/h a ±- 



A. Isotropic case 

In this subsection we briefly summarize the results ob- 
tained for the isotropic case in [10j. Numerical solutions 
of Eq. ^ derived for the isotropic case (A e s = 0) with 
various initial conditions are plotted in Fig.Q] (for a set 
of parameters given in the figure caption). Two fixed 
points, a repulsive (circle) one and an attractive one ap- 
pear in the phase diagram. Both can be determined by 
the analytical solution of the algebric fixed point equa- 
tion derived from ([7]) in case of vanishing anisotropy i.e. 
for A c ff = 0. The solution for the attractive fixed point 
is (see Eq. (25) in [13), 



UxO 

UyQ 



2o- 



2ua 



1-u 



(10) 



(in case of the repulsive fixed point u x0 and u z0 have 
been multiplied by -1). The attractive fixed point of the 
LLG equation in the rotating frame corresponds to the 








FIG. 1: Phase portrait in the rotating frame obtained by 
solving the LLG equation ((9J for the isotropic case (A e ff = 0) 
with the parameters cw = 0.1, w = 0.01, uil = 0.2. The circle 
indicates the repulsive fixed point at <f> = 3.12, 9 — 1.53. The 
attractive fixed point is at = 0.02, 6 = 1.61. 



stable steady state solution obtained in the laboratory 
frame (Eq. (24) in [13) 

M x (t) = u x o cos(ujt) — UyQ sin(wi), 
M y (t) = u x o sin(wi) + u y o cos(wt), 
M z (t) = u z0 . (11) 

The steady state solution enables us to calculate the en- 
ergy loss for a single particle. The energy dissipated in a 
single cycle can be calculated as (based on Eq. (fTTjl ) 



fi m s 



dt H- 



dM 
~~dT 



= ^2-KmsH (-u y o) , (12) 



(see also Eq. (32) in [10() which has the form in the low- 
frequency limit, uj <C ajv, 



E = 2ir/iQmsH 



a N uj 2 L uj 3 



a 



N 



(13) 



Let us note that in Ref. [T(| the energy loss per cycle of 
isotropic nanoparticles obtained by oscillating and rotat- 
ing external fields have been analyzed in the absence of 
thermal effects, see e.g. Fig. 2 in It was shown that 
in the low-frequency limit the energy loss per cycle was 
found to be larger in the linearly polarized case. In or- 
der to consider the role of thermal fluctuations in case of 
isotropic samples let us compare the findings of [f| (where 
thermal effects were included) to [lOj. Dashed lines on 
Fig. 6 of [H correspond to the limit T — > and agree to 
the findings plotted in Fig. 2a of [l(| qualitatively. The 
important result is that in the limit of low frequency, the 
linearly polarized field was found to produce more heat 
power both for T = and for T =/= 0. 
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B. Anisotropic case 

In order to study the role of anisotropy let us follow 
the strategy applied in the isotropic case. Namely, we 
calculate the loss energy per cycle by determining the 
attractive fixed point solution of ((J) for non-vanishing 
anisotropy (A e ff ^ 0) which reads 

= cuuy + a N u 2 + a N u 2 z - lu l \ cS u v u z - aNX c su x u 2 z , 
= —uju x 
= uj l u v 



- u L u z - a N u x u y + uj L \ cfi u x u z - a N X QS u y ul, 
a N u x u z + a N X cS (l - u 2 z )u z , (14) 



which can be reduced for the following equation for u z , 



ra N 



1 2> 



A 



efi 



1 



(15) 








with r = w/(cj| + ot 2 N ). Let us note that Eq. (Ti"5"j) is 
identical to Eq. (3.8) of [l2[ (with zero static field H = 
and with p = 1) if the following identifications are used: 
h = 1/A ff, k = rwi/Aeff and A = a at /uil- The switching 
of the nanoparticle magnetic moments and the dynamical 
effects under the action of rotating field have been studied 
in [l| with great details and further considered in [13, 
0j|, however energy losses were not calculated which is 
the goal of the present work. 

Eq. (fTS"]) can be solved analytically which is used to 
calculate all the Cartesian coordinates of the attractive 
fixed point (or fixed points) in the rotating frame. If 
u y o (y-axis component of the attractive fixed point) is 
known the energy loss per cycle can be evaluated simi- 
larly to the isotropic case. For biomedical applications 
such as hyperthermia the the low-frequency limit is rel- 
evant. Therefore, let us first consider the low-frequency 
oj <C o>n and small anisotropy A c ff <C 1 limit where one 
finds a single attractive fixed point with 



lyO 



a N uilu] 3 



a 



N 



4r) 3 



(l + 2A off ). (16) 



Inserting (|16[) into the expression of the energy loss per 
cycle (|12p one finds, 



a 



N 



4) 3 



(1 + 2A eff ) 



(17) 



which shows that in the low- frequency limit, the small 
anisotropy does not modify the energy loss per cycle ob- 
tained for the isotropic case. At higher frequencies in 
case of small anisotropy a decrease is observed in the en- 
ergy dissipated in a single cycle compared to the isotropic 
case. This analytic result is supported by the numerical 
integration of Eq. ([9]) , (low- frequency, small anisotropy) , 
see Fig. [21 which is very similar to the isotropic case, see 
Fig.[U Figures Fig. [T] and Fig. [5] are similar but not iden- 
tical. Trajectories of the two figures, far from the fixed 



FIG. 2: Phase portrait in the rotating frame obtained by 
solving the LLG equation © in the limit of low-frequency and 
small anisotropy for the parameters ajy = 0.1, w = 0.01, lol = 
0.2 and A e ff = 0.05. It is similar but not identical to Fig.[T] 



points differ from each other (but the deviation is small) . 
However, for low- frequency and small anisotropy, the at- 
tractive and the repulsive fixed points of Fig.[T]and Fig. [2] 
coincide. The stability analysis done in [l2| also confirms 
the existence of a single attractive fixed point in this 
regime of the parameter space. The difference between 
the positions of the attractive fixed point obtained for 
the isotropic and anisotropic cases is more recognizable 
at higher angular frequencies. For example, according 
to the approximate expression (TlTJ)) . the attractive fixed 
point of Fig.[3]is at u v q = —0.0956, numerical results give 
Uyo = —0.0957 and the corresponding isotropic case (see, 
Eq. (Ti"7j]0 gives u y o — —0.0961. Thus the energy loss per 
cycle (which is related to — u v q) is decreased for the set of 
parameters used in Fig. [3] compared to the isotropic case 
(for the same w, ujl and ajv)- 

Let us consider the low-frequency but large-anisotropy 
limit. In this case one finds two attractive fixed points, 
see Fig.|U The attractive fixed points are situated above 
(up) and below (down) the equator. In the strong 
anisotropy limit their ^-components are the same and 
their ^-components are symmetric to the equator, thus 
their Cartesian components are u"q 



XO J l U 20l 

u down xherefore, both attractive fixed 
points have the same y-components in the rotating frame 
which reads as 



|<g wn |and^ 



UyO 



1 



a 



(18) 



N 'Yff 



Inserting (|18[) into the expression of energy loss per cycle 
(Il2l one finds, 



E « 2TTfj,o m sH 



a 



N 



(19) 
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FIG. 3: Phase portrait in the rotating frame obtained by 
solving the LLG equation <(9j (medium-frequency and small 
anisotropy) for the parameters ajv = 0.1, ui = 0.05, ljl ~ 0.2 
and Aeft = 0.05. According to numerical results, the attractive 
fixed point is at cj> = 0.097995, 6 = 1.77908 which gives u y0 = 
-0.0957. 
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FIG. 4: Phase portrait in the rotating frame obtained by 
solving the LLG equation (|9]) in the limit of low-frequency and 
large anisotropy for the parameters cln = 0.1, uj — 0.01, ujl = 
0.2 and A c g = 2.5. Two attractive fixed points appear in the 
figure. The repulsive fixed point and the saddle point are 
indicated by a circle and a solid circle, respectively. 



which vanishes for A c ff — > oo. In the limit of extreme 
large anisotropy there is no room for energy dissipation 
since the magnetization is aligned to the easy axis inde- 
pendently of the applied field. Indeed, for X e g — > oo the 
xy-plane components tend to zero u x q — ¥ 0, u v q — ¥ and 
u z — > ±1. In the (4>-9) plane the attractive fixed points 



should tend to the "poles", i,e, = and 6 = 0, w. If 
the anisotropy is decreased, they "move away" from the 
poles and tend to the equator, see Fig. [3] A critical value 
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FIG. 5: Phase portrait in the rotating frame obtained 
by solving the LLG equation <(9j (low- frequency, large 
anisotropy) for the parameters qjv = 0.1, ui — 0.01, lol = 0.2 
and Aoff = 1.5. The attractive fixed points are closer to the 
equator as compared to Fig.|4]where the anisotropy was larger. 

for the anisotropy parameter can be identified where one 
of the attractive fixed point (the one which corresponds 
to small 8, i.e. which lies above the equator ©) van- 
ishes. The other attractive fixed point remains always 
below the equator. The phase diagram obtained in the 
low-frequency limit, slightly below the critical value of 
anisotropy is plotted in Fig.|BJ 

Let us consider the loss energy as a function of the 
anisotropy parameter. In Fig. [7] the loss energy is plotted 
versus A e ff (for a set of parameters given in the figure) 
and it shows that in case of a rotating field the loss energy 
is obtained to be a monotonic function of the anisotropy 
parameter. Thus, according to our results, the anisotropy 
(where the easy axis is perpendicular to the rotating ex- 
ternal field) cannot be used to increase the heating ef- 
ficiency of magnetic nanoparticles in the low-frequency 
limit. 

For the sake of completeness let us consider the high- 
frequency limit, although it is out of the scope of the 
present work (irrelevant in case of cancer therapy by hy- 
perthermia), hence, we do not study this in detail. The 
phase diagram obtained in the high-frequency and large 
anisotropy limit is shown in Fig. [8] and can be compared 
to the one obtained for low frequencies (with the same 
anisotropy), see Fig.|U The similarity between the high 
and low frequency cases is that two attractive fixed points 
appear. However, in the high-frequency case they have 
different coordinates in the (<fi-6) plane. Thus, the energy 
losses correspond to the attractive fixed points differ from 
each other. 
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FIG. 6: Phase portrait in the rotating frame obtained by 
solving the LLG equation ([9]) in the limit of low-frequency, 
slightly below the critical value of anisotropy. The parameters 
are ajv = 0.1, us = 0.01, uil = 0.2 and A e a = 1.175. There is 
only a single attractive fixed point below the equator (large 9). 
The other attractive fixed point (above the equator, small 8) 
of the large anisotropic case and also the saddle point vanish 
in this regime of the parameter space. 



FIG. 7: The energy loss is plottes versus the anisotropy 
A c ff. There are two scaling regimes, the one is at small, the 
other is at large anisotropy which are separated by the critical 
value Acr. The full line corresponds to energy loss obtained at 
the stable fixed point below the equator and the dashed line 
represents energy loss at the fixed point above the equator 
which is stable only for A c ff > A cr . 



In summary we conclude that the uniaxial anisotropy 
(where the easy axis is perpendicular to the rotating ex- 
ternal field) either does not modify the energy loss per 
cycle (in case of small anisotropy) or the energy dissi- 
pated is decreased as compared to the isotropic case. 



IV. LINEARLY POLARIZED APPLIED FIELD 
IN THE LIMIT OF LARGE ANISOTROPY 

One of the goal of this paper is to investigate the role of 
anisotropy in the possible enhancement of heat efficiency 
of magnetic nanoparticles driven by a rotating magnetic 
field. On the one hand, in the previous section it was ob- 
tained that in case of circularly polarized applied field if 
the uniaxial anisotropy (perpendicular to the rotating ex- 
ternal field) has been taken into account, the energy loss 
per cycle either remains unchanged (small anisotropy) 
or decreased (large anisotropy). On the other hand, in 
Ref. [10( it was shown that for isotropic nanoparticles 
the linearly polarized external field provides us a larger 
heat power in the limit of low frequency. Furthermore, 
the latter statement was shown [5j to be reliable in the 
presence of thermal effect, too. Thus, the circularly po- 
larized applied field cannot be used to achieve a better 
heating efficiency by nanoparticle systems (neither for 
isotropic, nor for anisotropic nanoparticles with uniaxial 
anisotropy) if the effect of thermal fluctuations is negli- 
gible. The role of thermal effects is discussed in IVl where 
it is argued that the possible modification of the above 



finding by thermal fluctuations can only be expected in 
case of moderate anisotropy. Thus, the results of the 
present work indicates that the heating efficiency cannot 
be increased by the rotating field at least for small and 
very large anisotropy in the limit of low frequencies (in- 
dependently whether thermal effect are included or not). 
This finding does not require any further analysis of the 
linearly polarized applied field. 

However, the study of energy losses in case of the lin- 
early polarized applied field for anisotropic nanoparticles 
enables us to consider whether the large anisotropy can 
be used to explain the experimental results of [7j . Let us 
note the formation of chains by nanoparticles is a known 
feature of ferrofluids and the chains of particles represents 
a very large anisotropy [2l[ . Therefore, in this section we 
discuss the solution of the LLG equations ([1]) obtained 
for a single-domain magnetic nanoparticle under linearly 
polarized, i.e. alternating applied (external) magnetic 
field in the limit of large anisotropy. The applied field 
is assumed to oscillate along the x-axis with an angular 
frequency us 

H cxt = ^ (cos(wt),0,0), (20) 

7 

where u>l = 7'|H| is the Larmor frequency. Similarly to 
the circularly polarized case, here, we consider particles 
with uniaxial anisotropy. The easy axis of the magneti- 
zation is chosen to be the x-axis (similar results can be 
obtained if it is chosen to be perpendicular to the x-axis) 

H aniso = ^ (A cff M x , 0,0), (21) 
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FIG. 8: Phase portrait in the rotating frame obtained by 
solving the LLG equation <(9j in the limit of high-frequency 
and large anisotropy for the parameters oin = 0.1, a; = 
0.15, ojl = 0.2 and A c ff = 2.5. There are two attractive fixed 
points in the figure similarly to the low-frequency case, Fig.[4] 
Here the attractive fixed points have different y-components, 
thus the corresponding energy losses also differ from each 
other. 



where M x is the x-component of the magnetization vec- 
tor and the parameter A e ff describes the strength of the 
anisotropy. Note that in case of alternating applied field, 
it is convenient to study the original LLG equation (Q} 
instead of the rotated one ([7]). The LLG equation for the 
Cartesian coordinates of the magnetization reads as 

^ = a N {\ cS M x + cos(W))(l - Ml), 
at 

— * = -{ui L M, + a N M x M y )(cos(ijt) + \ e gM x ), 
dM z 

= (cj L M y + a N M x M z ){cos{u)t) + X eS M x ), (22) 

which has in general no analytic solution. In the limit of 
extremely large anisotropy, cos(cjt) <C A e ffM x , however, 
the time-dependence of M x can be determined as 

Mx {t) = ± - eMaN\ cS t) _ (23) 
v /exp(2aArA off i) - C 

with C = 1 - 1/M% . If t -> oo the solution ([23} tends 
to ±1 and the energy loss per cycle vanishes. 

We conclude that in case of a very large anisotropy 
the heat power of a magnetic nanoparticle driven by a 
linearly polarized applied field vanishes. The same was 
observed in case of the circularly polarized external field. 
Thus, it is a natural requirement to obtain a comparable 
heat power given by the linearly and the circularly polar- 
ized applied fields if the anisotropy is large enough which 
can explain the experimental results of 7] . Indeed, if the 



fcrrofluid was not prepared appropriately, the nanopar- 
ticles can form chains and consequently, the anisotropy 
could become large and the energy loss tends to zero 
rapidly. 



V. THERMAL EFFECTS 

In this work we studied the influence of anisotropy on 
the energy losses in the framework of the determinis- 
tic LLG equation in the absence of thermal fluctuations 
for rotating applied field. In this subsection we discuss 
briefly how thermal effects (see e.g. |22H24| ) can pos- 
sibly modify the results obtained by considering purely 
dynamical effects based on the LLG equation. Let us 
follow Refs. [TH, [TJ, [l5j]. In case of a rotating field, the 
steady state solutions of the dynamical problem are inde- 
pendent of the initial conditions of the individual parti- 
cles. Therefore, the average magnetization can be easily 
determined by these (stable) steady states. 

In case of two steady state magnetizations (up and 
down states), due to thermal effects, a nonzero probabil- 
ity of a transition from one stable state to another ap- 
pears. This type of relaxation mechanism is missing in 
the present work due to the lack of thermal effects. We 
showed that large anisotropy is needed in order to have 
more then one stable state. It was also shown that if the 
anisotropy is not large enough only a single fixed point 
appear in the phase portrait, thus one has to consider 
only a single stable state. In this case the modification 
caused by thermal effects is less important. Thus, ther- 
mal effects can only modify the determination of energy 
losses (in case of a rotating applied field, for low frequen- 
cies relevant to hyperthermia) if the anisotropy is large 
enough but not too large (otherwise the barrier between 
one state to an other becomes too large). 

Furthermore, let us compare @ and in order to 
consider the role of thermal effects in case of isotropic 
samples. In these articles energy losses were investigated 
both for alternating and for rotating applied field. Ther- 
mal effects were included in [5j while these are absent in 
10]. For example, the dashed lines on Fig. 6 of Q corre- 
spond to the limit T — > and agree to the findings plotted 
in Fig. 2a of [l(| qualitatively. Let us pay the attention 
of the reader to the logarithmic scale used in Fig. 2a of 
[Icj . The modification caused by thermal effects (solid 
lines in Fig. 6) are less important for the rotating case 
but very significant for the alternating case in the limit 
of low frequencies. Nevertheless, for low frequencies, the 
linearly polarized field was found to produce more heat 
power both for T — and for 

According to Refs. [IH, [TJ, [l5[ for small and for very 
large anisotropy the thermal effects are less important. 
It was shown in Ref. § that for isotropic samples even 
for T ^= the energy absorption per cycle for a nanopar- 
ticle was larger in case of a linearly polarized applied 
field. Thus, if anisotropy (at least if it is small or very 
large) does not increase the energy losses obtained in 
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case of a rotating external field (which is one of the find- 
ing of the present work) , a better heating efficiency can 
be achieved for isotropic nanoparticles using linearly po- 
larized field. This indicates that possible modification 
caused by thermal effects can only be expected in case of 
moderate anisotropy. Therefore, the study of the present 
work is relevant for applications in hyperthermia at least 
for small and very large anisotropy. 

VI. SUMMARY 

The nonlinear dynamics of magnetization and the loss 
energy of a single magnetic nanoparticle with uniaxial 
anisotropy have been considered under circularly polar- 
ized applied field in the absence of thermal fluctuations. 
The easy axis of magnetization has been chosen to be per- 
pendicular to the rotating applied field. We solved the 
deterministic Landau-Lifshitz-Gilbert equation in order 
to determine the loss energy per cycle in case of the ro- 
tating applied field and the findings were compared to 
that of obtained for the isotropic case in [To|. Compari- 
son between the linearly and circularly polarized applied 
field has also been performed and the results were ana- 
lyzed in terms of the experimental data Q ■ 



Our goal was twofold: (i) to study whether the 
anisotropy can be used to achieve a better heating ef- 
ficiency in case of rotating external field, (ii) to use the 
anisotropy to resolve the discrepancy between theory (lpj 
and experiment [7j • We showed that for circularly polar- 
ized field, the loss energy per cycle is decreased by the 
anisotropy compared to the isotropic case. Thus in the 
low-frequency limit, more heat power can be achieved by 
alternating applied field for isotropic nanoparticles, at 
least the rotating applied field produces lower energy ab- 
sorption for small and very large anisotropy. It was also 
shown that in the limit of extremely large anisotropy, 
experimental results of [7J can be explained. The possi- 
ble role of thermal fluctuations discussed here indicates 
the necessity of the extension of the present study for 
the case of moderate anisotropy when thermal effects are 
taken into account appropriately. 
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